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ABSTRACT 

The 77 p^p^ An reaction shows a broad peak at 1.5 GeV in the 
{J^ , Jz) = (2^,2) channel which has no counterpart in the p^p^ channel. 
This "resonance" is considered as a candidate for a qqqq state in the "s- 
channel". We show, however, that it can also be explained by potential 
scattering of p°p° via the a- exchange in the "t-channel". 



1. Introduction 

After many successes of quark models to describe a single meson or 
baryon, their predictive power in the two-hadron sector remains question- 
able. A fair description of the nucleon-nucleon interaction |1]] was accompa- 
nied with the prediction of a rich dibaryon spectrum which has never been 
observed; it is not yet explained whether the predictions are wrong or the lev- 
els just too broad. Therefore it is of great interest to see how different quark 
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models perform in the two-meson sector, as compared to effective mesonic 
models. 

Much experimental and theoretical interest has concentrated in the pp 
system. Firstly, there is a remarkable experimental "puzzle". The 77 — > pp 
reaction is characterized by the following features 0, |^: a large and broad 
peak near the nominal threshold in the p°p° (2"*", 2) channel and a much 
smaller cross section ("noise") in other channels as well as in the p'^ p~ 
channel. 

Secondly, the pp system consists of two unstable (shortlived) hadrons. 
Therefore, the shape of the "resonance" is unusual and it challenges an ap- 
propriate quantum mechanical treatment. Also, such shape might be selec- 
tive regarding competing models. More experimental data would be very 
welcome in order to sharpen the details of this shape. 

A good model should explain all these features simultaneously. Since 
so far the quark models |^ covered more features simultaneously than 
the effective mesonic models P, ^, the opinion prevails that explicit quark 
degrees of freedom (a qqqq molecule) are essential for this process. In the 
present paper we challenge this view. We show that the considered process 
can also be explained by potential scattering. 

Of course, also in potential scattering, the intermediate system consists 
of two quarks and antiquarks, each qq pair being correlated like in a free p 
meson. This is just one important configuration, which also extends into the 
asymptotic region. Other intermediate qqqq configurations (such as those 
suggested in 0, |]) may or may not be important. They can easily lead to a 
resonance peak but this does not prove that they really dominate (that they 
are at the right energy and strongly coupled to the asymptotic wavefunction) 
unless no other mechanism exists. 

The peak near the nominal threshold and its strong spin dependence are 
reminiscent of the p-n scattering at low energy. There the S=0 channel has 
an exceptionally high cross section while in the S=l channel it is much lower. 
The reason is that in the S=0 state the potential well contains almost exactly 
1/4 wavelength of the relative wavefunction while in the S=l state there is 
a little more than A/4 in the well because the potential is slightly stronger. 
We hypothesize that a similar trick of nature is played also in the case of the 
p — p scattering. 

We assume that both photons are converted in p via vector dominance, 
and both p then interact by a scalar isoscalar potential which is just strong 
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enough to (almost) bind them for one spin orientation and fails to do so for 
the other spin orientations. In addition, an isoscalar potential has the desired 
feature that it cannot exchange charges and does not lead to the p^p^ final 
state. 



2. A nonrelativistic potential model 

In order to illustrate this idea, we choose a Yukawa-type a-exchange po- 
tential V{r) = Voexp(— mo-r)/(— mo-r). The full process is shown in Fig. 1 
and its Lorentz invariant amplitude can be factorized in the following way: 

M = M,,,, /pi. . (1) 

We treat the 7 ^ p conversion and p ^ tttt decay perturbatively, but 
the pp pp scattering nonperturbatively by solving the nonrelativistic 
Schrodinger equation. In the calculation of the matrix element for the process 
7i 72 — Pi P2 we account for inelastic scattering of p-mesons: 

M^^pp = J ipyj{f^i,f^2) Vjp i>ppirpi,fp2) df^idf^2dfpidfpi (2) 

where V^p is the contact 7p potential V^p = (f^p/M^) S^{f^i—fpi)6^{f^2—'f^p2), 
ijj^^ is the function of incident photons (plain wave, normalized in a box of 
volume V) and ippp is the wave function of the scattered p mesons, calculated 
by solving the Schrodinger equation. The amplitude for j —>■ p conversion 
f\p is determined from the experimental value for the partial width of the 
process p e^e~ [|] and is equal to f^p/m'^ = (2.01 * 137)"^. In eq. (||) the 
photon selects the correct momentum component (Fourier component) from 
the calculated scattering wave function of the two mesons. 

We account for the shortlived nature of p mesons (F = 152 MeV) while 
writing the propagators of the scattered mesons Ppi and Pp2: 



We reduce the Lorentz invariant four-body phase space (i$4 P to two- 
body phase space of the two p mesons by integrating over all pion momenta 
and assuming that the amplitude is not sensitive to them. Neglectig the 
contributions of order F^/m^ we get the cross-section 



3 



^ ^pir mp2T 

(m^]^ — m^)2 + ^pF^ (mp2 — w^p)^ + WT-pF^ ' 

where mpi and are the effective masses of the scattered p-mesons = 
Wp^—p^. This formula (used frequently before) has an intuitive interpretation 
that the two "final" p mesons have arbitrary masses with a Breit-Wigner 
probability distribution around the nominal mass rrip = 770 MeV. 

Since the scattering phase shift is very sensitive to potential parameters 
as well as on p masses we do calculate it separately for each combination 
of the two p masses and add up the contributions in the double integral 
numerically, with the appropriate Breit-Wigner weights. 

The calculated cross sections a{W) are very sensitive to potential strenght 
Vq and weakly sensitive to potential range 1/mo- (in this nonrelativistic model, 
rua = 550 MeV is taken). Results for different Vq are presented in Fig. 2. 

For a narrow range of the model parameter Vq, the results qualitatively 
agree with the experimental data, as long as we are considering the height of 
the peak in the (2"*", 2) channnel. The position of the peak comes, however, 
too high above the nominal threshold. Similar calculation using relativistic 
Klein-Gordon equation for the p-meson scattering on fixed potential shows, 
that the tendency of relativistic effect is to pull the peak down to lower 
energies. It will be, however, shown in Sect. 3 that a more severe calcula- 
tion including relativistic effects and gauge invariant coupling can cure this 
deficiency of the toy model. 

In order to get a very low cross section in the (2"'', 0) channnel one has to 
assume a strong spin dependence of the cr-exchange potential. Such strong 
spin dependence is arbitrary in the present nonrelativistic model, but it will in 
fact appear in the relativistic model in Sect. 3 as a consequence of momentum 
dependence of the gauge invariant coupling. 

3. A relativistic model with gauge invariant ppa coupling 

The Lagrangian, from which we obtain the potential, respects vector dom- 
inance and is gauge invariant: 
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PC-?F,,{B)F^'^{B)] 



(5) 



where F^y{^) 



d $ 



dy^^. After diagonalyzing the quadratic term in the 



Lagrangian, we may want to have a vanishing (T77 vertex, as approximately 
required by the low-energy 77 — > vrvr scattering and confirmed in some phe- 
nomenological models as the NJL model with a sharp cutoff 0. This fixes 
j3 = 1 — a. Such an assumption is however not needed in our calculation, 
since our result in the leading order of e/g does not depend on the parameter 
f3. 

The unitary amplitude is calculated with the Bethe-Salpeter equation 
with the kernel given by Eq. (j^). The equation is simplified by the ladder 
approximation (only a exchange in the t and u channels) and by taking into 
account only the leading contributions in powers of e/g. The Blankenbecler- 



Sugar reduction ||T^ is used (with no retardation, i.e. with equal off-shellness 
of the final and intermediate p) and the partial wave expansion ||TT 
performed: 
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dk 



J K, — q^ — It 
^A'iA'2A"iA"2 

(g',^)/A",A".(A:)Mif„ 

iA"2AiA2 (A;,g), (6) 



/a"iA"2 



(47r) Vfc2 + 
where tv?{x) = W'^/A — and the factor 



aA"i(m {k))ax",,{Tn (k)) 




g2 < and A = 
g2 = and A = 
otherwise 



(7) 



(8) 
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takes into account that the longitudinal polarization vector eQ{q) becomes 
timelike (lightlike) for < (=0). 

In eq. the initial particles with polarizations Ai and A2 have the 
linear momentum q and are on-shell {pf = p\ = W^/A — q^ = 0). The final 
(intermediate) p mesons with polarizations A^, A2 (A'/, A2) can be equally 
off-shell with pj = pj = - q'^ {W^/A - k"^). If the total energy W is 

larger than the nominal threshold 2m, then = — > and eq. 

(H) has complex solutions as required by unitarity. 

Taking special linear combinations of amplitudes with definite polariza- 
tions one gets amphtudes with definite parity. Each channel gets decoupled 
from the others and can be represented by a system of one-dimensional inte- 
gral equations. The system for the (2+,2) channel is explicitly given in Ap- 
pendix A. These equations were then solved with the matrix inversion method 
[|T^ , taking 24 points for the magnitude of the relative three- momentum from 
to the cutoff A. It was also checked, that the results change slightly, if one 
takes only 16 grid points (the differences get however larger with energy, that 
is why we limited our calculation up to W = 1.8 GeV). 

The cross section was then calculated by weighing the amplitude squared 
with two Breit-Wigner factors (one for each final p). The approximation of 
equal off-shellness of the two p, with which the Bethe-Salpeter equation (|^) 
was derived and solved, is equivalent to the approximation that the amplitude 
depends only on the sum of the "masses" of the final p. For the total cross- 
section we get 

^TOT = / dq'p{q')-Y.i'2J+l)\Mi'!U,^,^{q',q)\', (9) 

where a first integration is already hidden in 

p{q') = K{^fj^ dp'[rh\p')+m\q')~rh\Q)]^^^ 

FBw{{rh{q') + p')\ml)FBw{{rh{q') -p')\ml) , (10) 

where K = [{e/g){l — a/2)]^/(87r^). The relativistic Breit-Wigner factor is 
defined as 
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Fswip', ml) - _ r^)^+ m^r^(p2) ^^^^ 

with the same p^-dependent p-meson decay width |jT^ as in the experimental 
analysis 0: 

r (v') = r f^2w P^-4m^ a ^H-Aml) 

with rup = 770 MeV, Tp{ml) = 152 MeV and = 140 MeV. 

The cross sections in separate channels are obviously only part of the total 
cross section (^). The explicit expression for the (2+, 2) channel is given in 
Appendix A. 

We take for the vector dominance factor e/g = from g^/{47r) = 2.3 [|I^]. 
The remaining parameter a were fitted to the experimental data in order 
to get the correct height of the cross section (the position of the peak and 
the relative contributions of different channels are however not dependent on 
a). We have checked that the fitted value for a gives a small enough (10~^) 
branching ratio for p° 77r°7r° decay, which is not observed experimen- 
tally. The mass in the two-body propagator after the Blankenbecler-Sugar 
reduction would be m = nip = 0.77 GeV for non-decaying p, but can be 
slightly modified in the more realistic case of broad p mesons. We used 
the value m = 0.692 GeV. Since we are interested in the energy range 1.4 
GeV< W < 1.8 GeV and since m < W/2 the obtained amplitudes are uni- 
tary. We checked that the cross sections behave similarly for W > 2m, if a 
different value for m is chosen. It came also out, that the results are rela- 
tively insensitive to the choice of the a mass, so that we will show only the 
cross sections for rricr = 0.5 GeV. 

The results for the interesting channels are shown in Fig. 3, for different 
choices of the free model paramaters A, and a together with the exper- 
imental values. We checked also that there is a combination of the three 
free parameters in a reasonable range, which reproduces quite well the ex- 
perimental results: if we fix one of them, the other two are rather precisely 
constrained by the data. 



4. Conclusion 
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We have shown that there exist plausible sets of parameters for poten- 
tial scattering (for the ppa coupling) such that the peak in the (2"'", 2) 
channel at 1.4 — 1.8 GeV is reproduced while in other channels the cross 
section remains small. The parameters are plausible if compared to the cr- 
exchange contribution of the N-N interaction (for example in Bonn potential 



[p!5|]). The agreement does not prove that the mechanism is just potential 
scattering, but it presents potential scattering as a candidate until proven 
otherwise. Of course, we neither prove that resonance scattering through a 
qqqq molecule is wrong, we only show that it is not conclusive since other 
mechanisms exist. A finer distinction will be needed. 

The succes of the potential scattering relies on the choice of model pa- 
rameters which "hold" about a quarter of the wavelength of the scattered 
p in the "potential well", a situation similar to the singlet n — p scattering 
where there is also almost exactly a quarter wavelength in the well yielding 
a high peak at low energy. The nonrelativistic toy model illustrates this 
point qualitatively; the position of the peak comes, however, too high above 
the nominal threshold. Relativistic effects and gauge invariant momentum 
dependent coupling cure this deficiency; they also give a strong spin depen- 
dence of the (T-exchange potential needed to keep the cross-section low in 
other channels. 

This work was supported by the Ministry of Science and Technology of 
Slovenia. 



APPENDIX A 

In this Appendix we will write down some useful formulas for the cal- 
culation of the amplitudes and cross sections. We will concentrate only on 
one specific channel, namely, the most important {J^, Jz) = (2^, 2) chan- 
nel. Similar equations can be (and were) defined also for the other possible 
channels. 

First one has to project the potential to states with definite total angu- 
lar momentum. This was done integrating the potentials for pp scattering 
obtained from the Lagrangian (^ over the scattering angle ||11||,||12||: 
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n^^kA,A.(fc', = 2^ £ d{cos9)Vx',,,,,,,,ik', k) . (13) 

Using these quantities one could solve for each J the Bethe-Salpeter equa- 
tion (^ for the 4 X 9 = 36 unknown functions (of q') M_[/y^^^^. However, 
due to parity and Bose symmetry, not all these functions are independent. 
In fact, they are related by 



^r^-A'.-A.-A. (14) 

(X = M or y) so that the matrices M{q', q), M{k, q) in V{q' , q) in (|^) have 
only 3 independent elements for J = and 12 for J = 2, while V{q' , k) has 
5 independent elements for J = and 25 for J = 2 (the difference comes 
from the fact that the matrix elements of V{q' , k) can have also longitudinally 
polarized particles in the ket, while only transversal polarizations are allowed 
for the other, since — = 0). 

All the elements Xy^y^^_^^^{q',q) do not have good parity. Using the 
equations (|T^ and taking special linear combinations we can decouple the 
original Bethe-Salpeter equation (^, obtaining lower dimensional matrix in- 
tegral equations, one for each (J^, J^). 

For example, we get for ( J^, J^) = (2+, 2) the following integral equation: 
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where we defined two averages 



-^a;a^AiA2 ^ 9(^A'iA^AiA2 + ^A^A'iAiA2) ' 



and where 



' ^ ' (4x)»VFW ^ ' 

and a+(a;) = +1, a^{x) = +1, 0,-1 for a; = +1, 0, —1 respectively. 
The cross section in this same channel is then 



1 da(^^'^^ (2) x|2^U>(2) . / x|2 



p{q') dq' 



5{\M['l^_{q',q)\' + \M[\4q',q)\' + 
2|Mil4g',g)|2 + i|Ma_(g',g)|2) . (20) 
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FIGURES 



Figure 1: Feynman diagram for the process 77 —>■ p^p^ innnx. 

Figure 2: The experimental cross sections for the (2"*", 2) channel (dots) 
and calculated cross sections for three different potential strengths Vq: Vq = 
1400 MeV gives the highest cross section, while Vq = 800 MeV gives the cross 
section with the same height as the experimental one. 

Figure 3: Cross sections for various channels in the relativistic potential 
model with scalar exchange in comparison with experimental data. 
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